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Abstract 

We discuss a phase diagram for a relativistic SU{2) x Us{l) lattice gauge theory, with 
emphasis on the formation of a parity-invariant chiral condensate, in the case when the Us{l) 
field is infinitely coupled, and the SU{2) field is moved away from infinite coupling by means 
of a strong-coupling expansion. We provide analytical arguments on the existence of (and 
partially derive) a critical line in coupling space, separating the phase of broken SU{2) symmetry 
from that where the symmetry is unbroken. We review uncoventional (Kosterlitz-Thouless 
type) superconducting properties of the model, upon coupling it to external electromagnetic 
potentials. We discuss the role of instantons of the unbroken subgroup U{1) G SU{2), in 
eventually destroying superconductivity under certain circumstances. The model may have 
applications to the theory of high-temperature superconductivity. In particular, we argue that 
in the regime of the couplings leading to the broken SU{2) phase, the model may provide an 
explanation on the appearance of a pseudo-gap phase, lying between the antiferromagnetic and 
the superconducting phases. In such a phase, a fermion mass gap appears in the theory, but 
there is no phase coherence, due to the Kosterlitz-Thouless mode of symmetry breaking. The 
absence of superconductivity in this phase is attributed to non-perturbative effects (instantons) 
of the gauge field U{1) G SU{2). 
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1 Introduction 



There has been a great deal of recent interest in the dynamical symmetry breaking patterns 
of three-dimensional quantum gauge field theories, both from the pure particle theory 
standpoint |I|, ^, and as a tool for describing models of High Tc Superconductors 

H, |[. The gauge theories studied in those works have been either QED3 and variants 
of it Bii,or SU{2)0U{1) ii. 

From the condensed-matter view point, which motivated our approach to the subject, 
the key suggestion which lead to non-abelian dynamical gauge symmetry structure for the 
doped antiferromagnet, was the slave- fermion spin-charge separation ansatz for physical 
electron operators at each lattice site i 0: 



C[ I. ^' \ -tpl tpl l.\ Z2 Zi 



where Xa/3 are 'particle-hole' matrix- valued operators [0, c^, a = 1,2 are electron anihi- 
lation operators, the Grassmann variables ipi, i = 1,2 play the role of holon excitations, 
while the bosonic fields Zi,i = 1,2, represent magnon excitations [Q. The ansatz (|I]) has 
spin-electric-charge separation, since only the fields ipi carry electric charge. The holon 
fields ipa/3 may be viewed as substructures of the physical electron Xa/s 0, in close analogy 
to the 'quarks' of QCD. 

As argued in ref. the ansatz is characterised by the following local phase (gauge) 
symmetry structure: 

G = SU{2) X Us{l) X Ue{1) (2) 

The Ue{^) electromagnetic symmetry is due to the electric charge of the holons. In the 
absence of external electromagnetic potentials is a global symmetry (fermion number). It 
becomes local (gauged) after coupling to external electromagnetism. 

The local SU(2) sjTiimetry is discovered if one defines the transformation properties 
of the Zap and -^^^ fields to be given by left multiplication with the SU{2) matrices, and 
pertains to the spin degrees of freedom. 

The local Us{^) 'statistical' phase symmetry allows fractional statistics of the spin and 
charge excitations. This is an exclusive feature of the three dimensional geometry, and is 
similar in spirit to the bosonization technique of the spin-charge separation ansatz of ref. 
IU|. The presence of Us{^) allows the alternative possibility of representing the holes as 



slave bosons and the spin excitations as fermions. 

In the model of 0, this UsiX) is assumed strongly coupled, capable of holon pairing 
and (parity-preserving) mass-gap generation. The mass generation breaks chiral sym- 
metry, which can be defined in three-dimensional theories with even number of fermion 
species [|I|, as is the case of the model of ref. 0. However, as discusssed in |]11], ||, |^, this 
mass gap is not accompanied by any phase-coherence, given that the symmetry breaking 



is realized in the Kosterlitz-Thouless mode 1112 . 
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At this stage we would like to make an important comment, concerning the relativistic 
nature of the effective model discussed in From a condensed-matter view point, such 
relativistic systems would arise by a linearization about specific points on the fermi surface 
of the statistical system, such as nodes etc. In this respect it is worthy of mentioning that 
recent experimental tests |]13[ imply that the superconducting gap in the high-Tc cuprates 



is of (i-wave type, with lines of nodes on the fermi surface. It is the linearization about 
such nodes, in the flux phase for the Us{l) gauge field, that leads to Dirac spectrum for 
ho Ion excitations, with the fermi velocity of holes playing the role of the limiting ('light') 
velocity, as suggested above and in refs. 0, ^ . Then, as a result of the Kosterlitz-Thouless 
(KT) mechanism for superconductivity described in ^, a fermion gap opens at those 
nodes, which, due to the absence of a local-order parameter, respects the d-wave character 
of the superconducting state. 

The pertinent long-wavelength lattice gauge model, describing the low-energy dynamics 
around such d-wave nodes, assumes the form [Q: 



+ 



E tr \kziVijUijZj + h.c. (3) 

<ij> 

where J is the Heisenberg antiferromagnetic interaction, is a normalization constant, 
and Aij is a Hubbard-Stratonovich field that linearizes four-electron interaction terms in 
the original Hubbard model, and Uij,Vij are the link variables for the Us{l) and SU{2) 
groups respectively. The conventional lattice gauge theory form of the action (D is de- 
rived upon freezing the fluctuations of the Ajj field , and integrating out the (massive) 
magnon fields, z, in the path integral. This latter operation yields appropriate Maxwell 
kinetic terms for the link variables VJj, Uij, in a low-energy derivative expansion []5|, p!4|. 
On the lattice such kinetic terms are given by plaquette terms of the form H : 



E [f^sumi^ - TrV,) + (3u,ii){l - TrU,)] (4) 



where p denotes sum over plaquettes of the lattice, and (3us{i) = Pi: Psu{2) = P2 = 4/3i are 
the dimensionless (in units of the lattice spacing) inverse square couplings of the Us{l) 
and SU{2) groups, respectively 0. The above relation between the /3j's is due to the 
specific form of the z-dependent terms in (^, which results in the same induced couplings 
9su{2) ~ 9us{i)' Moreover, there is a non-trivial connection of the gauge group couplings 

toiri: 

K oc glu(2) = 9ls{i) ~ (5) 
with rj = 1 — 5, 5 being the doping concentration in the sample P, To cast the 
symmetry structure in a form that is familiar to particle physicists, one may change 
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representation of the SU{2) group, and instead of working with 2x2 matrices in (P, 
one may use a representation in which the fermionic matrices ipap are represented as 
two-component (Dirac) spinors in 'colour' space: 



= (V^i - V-l)^ , ^,^ = (V^2 ^l)^ , i = Lattice site (6) 

In this representation the two-component spinors \E' (|^) will act as Dirac spinors on 
the Lattice, and the 7-matrix (space-time) structure will be spanned by the irreducible 
2x2 Dirac representation. By assuming a background Us{l) field of flux tt per lattice 
plaquette 0], and considering quantum fluctuations around this background for the f/s(l) 
gauge field, one can show that there is a Dirac-like structure in the fermion spectrum |[I6|, 



T7[ ^ |T5[, which leads to a conventional Lattice gauge theory form for the effective low- 
energy Hamiltonian of the large-f/, doped Hubbard model 0. 

In the above context, a strongly coupled Us{l) group can dynamically generate a 
mass gap in the holon spectrum 0, which breaks the SU{2) local symmetry down to its 
Abelian subgroup U{1) generated by the matrix. From the view point of the statistical 
model @, the breaking of the SU{2) symmetry down to its Abelian subgroup may 
be interpreted as restricting the holon hopping effectively to a single sublattice, since the 
intrasublattice hopping is suppressed by the mass of the gauge bosons. In a low-energy 
effective theory of the massless degrees of freedom this reproduces the results of ref. |^ 



derived under a large-spin approximation for the antiferromagnet, 5* ^ 00, which is not 
necessary in the present approach. 

The Kosterlitz-Thouless (KT) nature of the Us{l) induced mass gap (absence of 



local order parameter), is a characteristic feature of gauge theories in 2 + 1 dimensions, as 
argued in |TT|, ^ . When applied to our non- Abelian model it leads to unconventional 
KT superconductivity, provided the gauge boson of the unbroken U{1) G SU{2) is mass- 
less. Due to the compactness of the U{1) gauge group, however, which is a distinctive 
feature of the non- Abelian gauge group nature of the spin-charge separation (|1|), there are 
non-perturbative effects (instantons), which are responsible for giving the gauge boson 
UiX) a small but finite mass |T^. This spoils superconductivity, leaving only a phase, 
characterised by pairing among the holons, without the existence of phase coherence. It 
is one of the points of this article to argue that such a phase may provide a possible expla- 
nation of the so-called 'pseudogap' phase of high-temperature superconductivity pO|, an 
intermediate non-superconducting phase, lying between the antiferromagnetic and d-wave 



superconducting |TH] phases. A preliminary discussion on this issue appeared in ref. [21 



At this point, we would like to mention that other authors have also used relativistic 



fermions to describe the underdoped or pseudo-gap phase of high-Tc materials [23, 23 



Their approaches, however, are different from ours: In ref. [^ , relativistic charge excita- 
tions are used as in our model [^ . Their relativistic nature is due to the adopted scenario 
that the fermi surface of the underdoped cuprates consists of four small pockets, centered 
around (±f , ±f ) in momentum space. However, the low-energy model used in that work. 
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and the nature of the gauge symmetries involved, are different from our modeL In the 



'nodal liquid' approach of ref. on the other hand, the relativistic Dirac-fermion ex- 
citations around the four nodes of the putative fermi surface in the underdoped situation 
are neutral, and, hence, from our point of view they correspond to spin degrees of freedom 
rather than holons. This leads to a different physical scenario for the pseudo-gap phase 
than the one discussed here and in ref. [ PT[] . 

In this article we shall discuss in some detail the phase structure of the SU(2) ® UsiX) 
gauge theory. Despite the above motivation from condensed-matter physics, the analysis 
and the techniques used will be those of Particle Physics, thereby making the results even 
applicable to Particle Physics applications of three-dimensional Gauge theories, such as 
Early Universe studies, or high-temperature field theories. In this respect we mention the 
work by Volovik , which pursues the analogy between the Physics of superfluid Helium 
and that of the Early Universe, in an attempt to suggest condensed-matter experiments 
that could shed light in the physics of an early stage of our Universe. We hope that our 
work in this article will serve the purpose of pointing out yet another condensed-matter 
example, that of high-temperature superconductors, which may be connected to Particle 
Physics. 

The structure of the article is as follows: in section 2 we review the basic symmety 
properties of the Lattice action, including a discussion on the issue of spontaneous and/or 
dynamical breaking of parity in the context of the applicability of the Vafa-Witten 
theorem on the lattice. In section 3 we derive part of the phase diagram of the SU (2) ® 
f/s(l) model, in the strong-coupling regime of the SU{2) gauge group. The analysis is 
a non-trivial application of standard Lattice strong-coupling expansions ^ to our 
model. In section 4 we review briefly the unconventional superconducting properties of 
the system upon coupling it to external electromagnetic fields. Emphasis is placed on 
the Kosterlitz-Thouless type of breaking of the electromagnetic symmetry, which is not 
accompanied by phase coherence. In section 5 we discuss briefly the role of instantons 
in destroying superconductivity, but maintaining pairing and fermion gap formation. In 
section 6 we discuss a possible application of the model to the physics of high-temperature 
superconductors, with emphasis on the abovementioned role of instantons in inducing a 
pseudo-gap phase. This is an exclusive feature of the non-Abelian model of [0. The 
possibility of tuning the doping concentration in the sample to reach supersymmetric 



points in coupling-constant space |]2S[, with interesting consequences, is also mentioned 
briefly. Moreover, the role of additional four-fermion interactions, which may dominate the 
superconducting phase, is pointed out. Conclusions and outlook are presented in section 
7. Some technical aspects of our approach, such as rules of strong-coupling expansion, and 
the evaluation of a Jacobian in the transition from fermionic Lattice variables to mesonic 
fields, are given in three Appendices. 
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2 Basic Symmetry Structure of the SU{2)^Us{l) the- 
ory 

2.1 The Lagrangian of the model and its symmetries 

The theory (||) corresponds, after integrating out the magnon degrees of freedom 2; to 
a (low-energy) lattice lagrangian given by p 



i.fi 

A - trUp) + E(l - trVp) (7) 
p p 

where t/j,^ = exp{i9i^fj_) represents the statistical Us{l) gauge field and V^,^ = exp{ia"'Ba) 
is the SU{2) gauge field. The quantity K = K\tij\, with \tij\ = t assumed small @]. The 
fermions are 2 component spinors in both Dirac (Greek) and colour (Latin) space, \E' = \E'°, 
and the generators of the SU{2) group are the 2x2 Pauli matrices, af^, i = 1,2, 3. The 
Dirac matrices can also be taken to have the Pauli matrix representation (we continue to 
write them as 7^^, to distinguish them from the SU{2) colour matrices). Here we have 
passed onto a three-dimensional Euclidean lattice formalism, in which is identified with 
For completeness we mention that the (naive) continuum lagrangian corresponding 
to (^ is given by: 

c = -l{F,.r - l{g,.r + ^1^,7,^ (8) 

where = — igia^ — ig20''^Ba^^, and F^^, Q^^ represent the field strengths for the 
Us{l), SU{2) gauge groups respectively. 

There are two sets of bilinears which transform as triplets under a SU{2) transformation: 

Ai = -i["^l^2 - '*'2^l] , A2 = i["^1^2 + ^^2^1] , ^3 = "^1^1 - "^2^2 ; 

J^3^ = ^la^^l - ^2C^M*2 (9) 

where = ^i"7o'^ and two SU{2) bilinear singlets given by, 

A4 =^1^1 + ^2^2, J^4,^=^icrf,'^i+^2cr^,^2 , 12 = 0,1,2 (10) 
In this approach one can define meson states ^: 

Mab,af3 = ^b,l3^aa (11) 
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where Latin letters a, b denote colour SU{2) indices, and Greek letters a, j3 denote Lorentz 
spinor indices. One can re-express the meson state M, which is a 2 x 2 matrix in both 
colour and Dirac space, in terms of the above bilinears p|; 

M= A3I.CT3 + Ail.ai + A2l.a2 + Ail.l + 

^3,/.7'^or3 + J'i,f.r-cri + T2,^l^(y2 + ^4,m7^1 (12) 

The first matrix written is in Dirac space, and the second is in ^SU{2) colour' space. 



CO 








Figure 1: Phase diagram for the SU(2) ® U{1) model. The solid line is the critical line 
which is determined in this work, separating the phases of broken SU{2) gauge symmetry 
from the phase where the symmetry is unbroken. Its precise shape is conjectural at present. 
Analytical and continuity arguments in this work determine the shape of the line in the 
neighborhood of {Pi, P2) = (0,0) and {Pi, P2) = {Pi, 00) only. This critical line also seems 
to characterise a solid state model, whose low- energy continuum limit is the gauge theory 
studied in this work. 

The interesting feature of the SU {2) xUs{l) model is that the parity- invariant conden- 
sate transforms as a SU{2) triplet (H), and, hence, once formed, it breaks SU{2) U{1) 
dynamically [0, p. The parity- violating condensate, one the other hand, is an SU{2) 
singlet. In continuum theories, the energetically preferable configuration in the absence 
of external sources is the parity- invariant condensate, according to the theorem of Vafa 
and Witten [^, ^ on the impossibility of spontaneous parity breaking in vector like the- 
ories, which we shall discuss in the next subsection. Thus, at least from naive continuous 



6 



considerations, one expects that Energetics favours the formation of parity invariant con- 
denates, and this was the main reason why parity violating condensates have been ignored, 
so far, in the existing hterature. As we shall argue in the next subsection, this feature is 
respected by the lattice model of [^. 

All these ideas can be incorporated into a rough phase diagram for a three-dimensional 
SU{2) X U{1) theory, proposed in [^. The diagram is depicted in fig. |I|. The couplings 
shown are inverse couplings, P2 = ^/o-qIi Pi = l/'^fi'i; where a is the lattice spacing. 

The top line /?2 = 00, /3i 7^ 0, corresponds to the QED^ case. For QED^ it is now 
generally accepted that there exists 0, ||, ^ a critical number of fermion flavours, below 
which there is dynamical formation of a chiral condensate and chiral symmetry breaking 
[jl], 0]. In the language of an effective theory, where the dimensionless coupling is taken 
to be the inverse of the number of fermion flavours , we can say that there is a critical 
coupling above which there is symmetry breaking. 

In figure |l|, the (inverse) critical coupling on the lattice is denoted as (31. The shaded 
area shows the weakly coupled SU{2) breaking, and the fact that we have no breaking 



on the /?2 = (ie infinitely coupled) SU (2) line, as discussed in and will be reviewed 
below, means -by continuity- that one can draw a tentative critical line separating the 
broken and unbroken phases. The issue of whether the point where the line hits the /32 = 
axis is at the origin or at a finite value of f]!^ is one which we shall resolve here. 

We note at this point that, in the context of our statistical model 0, there is the 
special relation (|^) among the (inverse) couplings of the SU{2) and Us{l) factors, namely 
(32 = 4/5i, which, as we have mentioned, orginates from the special structure of the magnon 
{CP^) degrees of freedom of the model. This special relation is interesting in that, when 
combined with the fact that the gauge couplings in the statistical model depend on the 
doping concentration of the superconducting system, implies the existence of extreme 
values for the doping concentration, above or below which the broken-S'[/(2) gapped phase 
is lost. As we shall argue in this work, the critical value /^l''** = 0, which implies that 
in the context of the present effective theory one cannot see a minimum coupling below 
which the SU{2) symmetry is restored. It is understood that in the condensed matter 
context such a minimum coupling, appropriate for the onset of antiferromagnetism, arises 
from the magnon CP^ sector. We shall discuss such issues in more detail in section 6. 

2.2 Parity and Fermions on the Lattice 
2.2.1 The Vafa-Witten Theorem in the Continuum 

Before embarking into a detailed analytical study of the phase structure of SU (2) (g) Us{l) 
theory, we would like first to devote some time on the important issue of parity symmetry 
for Lattice gauge models. 

As is well known, in continuum models, an important theorem, due to Vafa and 



Witten forbids the spontaneous breaking of parity symmetry in vector-like theories. 
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in the sense that the parity-violating condensate is not energetically preferable. Let us 
briefly review this, in the context of our three-dimensional gauge model ||^. We shall 
consider the Euclidian path integral for the two different mass terms, corresponding to 
the condensates ^3 (parity preserving), and A4 (parity violating). Let us start from the 
case where the < Aa, >^ condensate is formed. In this case, the relevant path integral 
reads: 

= y DAD'^D^exp J d^x[L[A] + + im)^] = 

J DA det[i^ + im] exp J d^xL[A] (13) 

where L[A] denotes the pure gauge part of the lagrangian. 

We can see that det[i p] is positive because given i = A\E' then i ^(72(72^1'*) = 
A(72cr2^*). Thus every eigenvalue is repeated twice and the determinant (the product of 
the e-values) is therefore real and positive. The gamma matrices are in Dirac space and 
the sigma matrices are in colour space. 

However, det[ip + im] is not real for the following reason: The eigenvalue eiquations, 
in the presence of the mass m read: 



(z^ + zm)^ = (A + im)^ , and (i^ + im) (72(72^*) = (A + zm) (720-2^*) (14) 

The two equations have the same eigenvalue, however on squaring each eigenvalue we 
get a complex number and therefore the determinant is complex. 

Let us now come to the case where the parity- invariant condensate is formed < ^3 >y^ 
0. In this case, the effective action reads: 

Za:, = J da det[ip + imaa] exp J d^xS[A] (15) 
Applying the same method again we now get 

{ip + iriKTs)'^ = (A + im)"^ , and {ip + zmcr3) (72(12^*) = (A - zm) (720-2^*) (16) 

Now we see that the eigenvalues come in complex conjugate pairs, and therefore the 
determinant is real and positive. 

Thus, the determinants in both cases have the same absolute value, but the deter- 
minant in the case of a parity- violating condensate has an extra phase factor making it 
complex. It is then straightforward to argue that the parity-violating case will not be 
energetically preferable [^. To this end, we note that the vacuum energy (in a box of 
volume V) is given by. 
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e-^^^ = J rf/i(A)e/'^'"^[^]rfet,- (17) 

where detj,j = 3, 4 denotes the result of the fermion determinant in the case where the ^3 
or condensates are formed, respectively, and dfi{A) is the measure for the gauge field 
integration and is positive; the same is true for the exponential e-f "^^^^I^J in the Euclidean 
formalism. By a generic result in complex integration calculus, then, the phase in det4 
can only make the integration smaller than that for det^, and therefore the vacuum energy 
associated with mass < ^4 > is larger than the vacuum energy for < ^3 >. So we can 
say that the energetically preferred mass term is the parity conserving one, and this is 
essentially the theorem of Vafa and Witten . Caution should be expressed in applying 
the theorem to the case of dynamical mass generation, due to the absence of bare fermion 
masses, which leads to the existence of fermion zero modes, that make the Dirac operator 



ill defined, in need of regularization. However, the rigorous analysis of ref. ||2^ deals with 
that case too. 

2.2.2 Wilson fermions and the breakdown of the Vafa-Wittem theorem on 
the lattice 



On the lattice, however, the issue is non trivial, and still unsettled. As argued recently pT 



although the Vafa-Wiiten theorem p5| may hold in the continuum limit, however, on the 
Lattice there may be terms (at least in an effective Lagrangian level), proportional to the 
lattice spacing a, which may violate explicitly the parity symmetry, thereby acting like 
external sources and hence spoiling basic assumptions of the Vafa-Wiiten theorem P^ . 
The issue of how the continuum limit is taken is therefore a tricky one, and currently 
there is a debate as to whether spontaneous breaking of parity occurs on the lattice 



32| . Although we shall not enter this debate, which concerns Wilsonian fermions on 
the lattice that we do not use here, however we consider it as useful to point out the 
difficulties associated with the parity symmetry, since it is a very important issue for the 
superconductivity mechanism of the models 0, ^. This will help the reader appreciate 
better how these problems are avoided in the specific lattice model of ref. [^], which we 
use for our purposes in this work. 

The main problem with the lattice formulations of the Vafa- Witten theorem, using 
Wilson fermions, is associated with the fact that, in the case of spontaneous breaking 
of parity, the Dirac operator has zero modes, as we shall discuss below, and thus needs 
regularization. Such a regularization is provided by adding appropriate sources (which 
may trigger parity breaking) in the effective action and then removing them. The presence 
of a source term violates the vector-like nature of the regularized theory, and in general 
the problem arises from commuting the limits of removing the source or sending the bare 
mass term to zero. 

Let us first review the situation in the case of four-dimensional gauge theories. The 
reduction to three-dimensional gauge theories with even number of fermionic species will 
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be straightforward, as we shall argue below. In the case of lattice regularization with 
Wilson fermions, the appropriate hermitean operator is not the Dirac operator but the 
overlap Hamiltonian 75W^(mo), where mo is a bare mass term needed for regularization 
of the Wilson-Dirac operator W. This operator is known to have fermionic zero modes. 
The latter lead to a non-zero spectral density of eigenvalues p(A,mo) around A = 0, in 
the limit of zero (bare) fermion mass. If p denotes the density of the Dirac operator —ip, 
then the following result holds : 



p{X,m) = .J"^^ p(VA^-m^) |A|>m 
V A"' — rri'' 

= |A| < m (18) 

As the fermion mass m tends to zero, the operator p(A, m) —>■ p(A) non uniformly. 

To trigger numerically spontaneous breaking of a symmetry one adds an appropriate 
symmetry-breaking source term and then removes it. In the case at hand, one should 
add a source term of the form ih'^'j^as'^ . As noted in ||32|, then, the parity- violating 
condensate, proposed to occur in Wilson fermions is proportional to p(0, mo) as the 
source term is removed, /i — > 0^: 

< i^75(T3^ >= T27rp(0, mo) (19) 

where mo is a bare mass term which should be removed, and for simplicity we assumed 
two fermion flavours. 



The debate in the current literature [|3T| , p2| concerns the ordering of the limits mo 



0, /i — s> in the case of Lattice theories with Wilson fermions. The presence of a non-zero 
physical mass mo renders the limit /i — >■ safe, the parity-violating condensate vanishes 
on both the lattice and the continuum formalisms, and, thus, there is no problem with the 



theorem of 25 







In three-dimensional lattice gauge theories, with an even number of fermionic species, 
as the models we are interested in, chiral and parity-symmetry breaking may be studied 
in full analogy with four-dimensional gauge theories, provided one works with a 4 x 4 
reducible Dirac representation generated by the following 7 matrices: 



In such a case, there are two matrices that anticommute with the set of the 7-Dirac 
matrices 

Chiral symmetry is then generated by 75, and is broken by the parity-invariant condensate, 
which in four-component notation for the spinors ^ reads: ^3 =< \E'\E' >. On the other 
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hand, the parity-violating fermion condensate is given by: A4 =< \E'A^ >, with: 



A = ^7375 = _\ j (22) 
It is, then, straightforward, following ||32|, to show that A4 obeys a relation of the form 



(p!9| ) for the three-dimensional case. 

We now remark that, in our condensed-matter inspired case (^, the fermions describ- 
ing holon nodal excitations in the d-wave state are Dirac Spinors (^ and not Wilson. 
According to the discussion in [Q, for such a case there are no consistency problems or 



ambiguities, as far as results in the continuum are concerned |2^. To show this, in our 
case, one should first note that, as we shall discuss in more detail in the next section, the 
effective potential for the meson fields M (|TT]) assumes the following generic form ^ 



Veff ~ Tr 5: (^AlnM, - J] P{M,{-^,)Mh,)j (23) 

where A is a numerical constant, P{x) denotes an appropriate polynomial in x, i is the 
lattice site, and Tr is taken over the (reducible) 4x4 Dirac indices. As mentioned, above, 
to study spontaneous parity breaking numerically, one should add to ([23| ) an appropriate 
source term: 

Vs = hM, (24) 
Following ref. P^, we assume the following form for the vacuum wave function of Mf. 



Mi = Ue'^^ = U {cose + iAsin^) (25) 

where A has been defined in (|2^) . 

From the specific M-dependence of P{x) in (^), one observes that the only dependence 
on 6 comes through the source term: 

Veff [h] ~ 4 (^hUcosO + A\nU - 3P{U^)) (26) 

By extremizing ( PB| ) with respect to 6, one obtain only the trivial minimum = 0, in 
constrast to the case of Wilson fermions, where the possibility for a non-trivial solution 
for 6 exists, due to the Wilson parity-breaking term . This solution implies that parity 
cannot br broken spontaneously, in agreement with the (continuum) theorem of Vafa and 



Witten 25 



Hence, for our purposes in this work from now on we assume that dynamical mass 
generation in our SU{2) x Us{l) model selects -due to energetics - the parity-invariant 
combination, which is accompanied by the presence of Goldstone Bosons due to the break- 
ing of Sf/(2) ^ f/(l) i|, 0,11. 
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3 Study of the Phase Diagram of the 3-d SU{2)0Us{l) 
Lattice Gauge Theory 

3.1 The Phase Diagram in the regime (32 = oo, /3i arbitrary 

In this regime of the phase diagram the system is equivalent to a strongly coupled QED^, 
with global SU{2) symmetry. This global symmetry acts like a flavour symmetry, and is 
represented by an index /. Because of this limiting procedure to QED^, the spinors will 
be kept two component (||), so that a smooth transition to the case of non-zero SU{2) 
coupling is guaranteed. The pertinent path integral is: 

Z = Y[J dUi,^d^(d¥i expi-S[^, ^, f/]), (27) 

i,fi 

where S denotes the action written as 



+(3iY.i^-trUp) . (28) 



Putting Pi = 0, we can do the Us (I) integral immediately 



i,/! f i,fj. 



nn^o^'(2v?/L) • (29) 



= :^Tr[M/7^M,+^7^] . (30) 

The quantity (m/ )°^ denotes the meson field (making the Dirac indices explicit) (\l'f )"(\E'f )'^. 
Iq^ denotes the zeroth order Modified Bessel function truncated to an order deter- 
mined by the number of the Grassmann (fermionic) degrees of freedom in the prob- 
lem 1^, In our case, due to the Dirac indices and 2 flavours of the lattice spinors 
one should retain terms in Iq^ up to O^y"^): 

/-(2V^) = l + 2/., + f + § + (31) 
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The path integral is flavour symmetric Z = Yl^ Z f and, hence, we may factor out and 
ignore this dependence. We wish to obtain a path integral for the meson field M/, which 
necessitates the inclusion of the Jacobian for the pertinent field transformation. This is 
calculated in Appendix 2, following [E6|. The result for the partition function reads: 



Zf = llj dMlexp[~Y.^ogdetMi + J2logl'J{2^^)] 

i i,fi 

= llJdMi explY: -I logdetM/Ml^^ + \ogll^{2^^)] . (32) 



The effective potential depends on the variable y defined in (|30[) . To determine its form, 
in terms of the condensate, one writes the VEV, < M/ >= m^W where 1^ is a unitary 
matrix. So DetM- ~ m? and yi^^ ~ K'^m? /2. 

From the discussion in the previous section, we know that the energetically preferable 
configuration is the parity-conserving one, in which half of the fermion 'flavours' acquire 
masses +m, and the rest acquire masses —m |2^, |l|]. Thus, the tree- level effective potential 
becomes, neglecting overall factors, 

■^■^3^ 2 16 72 3072 ' ^ ^ 

which is plotted in fig. ^ (with = 1) [|. 

This potential has stationary points and, from the arguments of (described in 
Appendix 2), this is sufficient to show that < M/ >^ 0. 

The incorporation of the gauge interactions will change the situation, and induce 
non-trivial dynamics which may result in a change in symmetry for some regions of the 
gauge coupling constants. From the previous result (^3]), and the discussion in section 
2, it becomes clear that for weak SU{2) and strong enough Us{l) the SU{2) symmetry 
is broken down to a U{1) subgroup 0. The non-trivial issue here is whether there exist 
critical (inverse) couplings f3^,i = 1,2, above which the symmetry is restored. According 
to earlier analyses, either in the continuum or on the lattice [jl], |30|, |3^, there appears to 



be a critical coupling /?f on the axis P2 = oo, Pi = free, above which dynamical mass 
generation due to the Us{l) group cannot take place. This is depicted in figure |^. 

In the large-A^ continuum theories Jll one usually identifies 8gfN = a, where a is 
kept fixed as — > 00, and plays the role of an effective ultraviolate cut-off. In our 

^The general case K ^ 1 amounts to adding an irrelevant iiT-dependent constant to the expression 
for the effective potential. This is an exclusive feature of the minimal gauge model (|^), because in 
that case one may absorb K in the normalization of the fermion fields. However, in the presence of 
additional non-miminal fermion interactions, e.g. four-fermion interactions, which, as we shall see, may 
characterize realistic models of doped antiferromagnets in their superconducting phases, the constant K 
can no longer be absorbed in a normalization of the fermion fields, and hence its magnitude acquires 
physical significance. We shall address such issues in section 6. For the purposes of this section, the 
minimal gauge model (|^) will suffice. 
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Figure 2: The effective potential in the case (32 = oo, (3i = 0, which coincides with the case 
of strongly coupled QED^ . The potential has stationary points, implying a non-vanishing 
chiral- symmetry breaking condensate. 

lattice action we may identify the inverse of our lattice spacing with a/8, in which 
case j3l ~ N^.., with A*";, the critical number of (four-component) fermion flavours, below 
which dynamical mass generation due to Us{^) occurs. To leading resummation 
Nc ~ 32/7r^; incorporating l/N"^ corrections shifts this number slighly higher. The issue on 
the existence of a critical number is still not quite settled, and proper lattice simulations 
are needed in this respect. For our qualitative purposes, however, the large-iV critical 
number will be sufficient. 

3.2 Strongly- Coupled SU{2) ® Us{l) regime 

In this subsection we commence our analysis of the effects of strongly-coupled SU{2) 
gauge interactions, by means of a small (32 expansion. 

3.2.1 The Effective Action for /^s = /?i = 0. 

First, we examine the model at the (limiting) point (32 = (3i = Q (ie the origin of the 
diagram of Fig. We absorb the paramerter 7^ in a redefinition of the fermion fields, 
because the action under consideration is only quadratic in ^ fields 0. In this case, one 

^This is because our lattice fermions are of Dirac type. In the Wilson fermion case, on the other 
hand, due to the Wilson and bare mass terms, that one is forced to add in the lattice action, the above 
normalization of the spinors by K cannot be done. 
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may integrate out first the Us{l) gauge field. The SU{2) action, then, is separable into 
an integral on each link on the lattice [p9| : 

/ rfy 1 1 + tr(AV)triAV^) + M'^VM^^')]" + [tr(AV)tr{AV^)f ^ [tr(AV)tr{AV^)r \ 

(34) 

where the variable A are defined as follows: 

A^ixT, = ^,{x + a)7^^'^(x), A^{x)t = ^fe(x)(-7;.)^"(x + a) . (35) 

The evaluation of these terms was first done by Samuel p^, whose formalism we follow 
here. The resulting partition function, now with all gauge fields integrated out, is: 



= n/ d'^id'^'^ n {l + + lMA,^^A^^.)f - ^tr[{Ai^^A,^^y] + 

~^[t^iA,tiA,ti)]'^tr[{Ai^ij^Ai^fj,y] + Y^tr[(ylj^^74j_^)^]^| . (36) 

Since we want Zq to be in the form Zq = e^"^^ we exponentiate the above polynomial still 
keeping terms up order 0{AA)'^. 



^0 = n/ dmidm^e^vY. \^tr{Ai^^Ai^^) + ^[tr{Ai^^Ai^^)f - ^tr[(A,,^Ai,^) 



+^^^(A,MA^,M)^r[(Z,,^A,,^)2] + -^[tr{\^.A, 



5760 



We want to rearrange the (\1', \E') which make up the A, A to get an effective action in terms 
of meson fields defined on site. The standard procedure for making an effective action in 
terms of meson fields is to rewrite each term in (p6D function of M'^"-"^ = ^'^'"^"'^ 
eg; 



= -trM'rtrM!U-l,T^{l,V' ■ (38) 
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tVcM means we have traced over the colour indices. However if we were to substitute in 
for M the expansion in terms of bihnears (p!^ ) we would find that the term containing 
^3 vanishes because of the traceless property of the Pauli matrix a^. Thus the above 
rearrangement of tr{AA) will not be of any use if we want to get information on the VEV 
of .43. We must thus look for an alternative rearrangement. 

This alternative arrangement can be demonstrated if we look at the term [tr(Aj^^Aj^^)]^, 
which is written out as; 



we can rearrarange the \1') which appear into pairs of meson states defined on each 
site. 



This time when we substitute in our bilinear expansion for Mf ° we do not lose the part 
depending on ^3. We are no longer tracing over a single Pauli matrix, but rather over 
tr{a^) = tr{l). This procedure will be justified if we find that there exists an energetically 
favorable non-zero VEV for ^3. 

Also we see that the terms in the action which have an odd number of (\E', pairs 
(ie. 0{AA) and 0{AA)^) vanish in this procedure since they cannot be arranged to solely 
produce pairs of mesons. There is a subtlety involved in calculating the 0{AA)^ term, 
which we describe in Appendix 1. 

Having followed this procedure for each term in the action the next step will be to 
substitute in for M, or more precisely for ^3, the VEV we are looking for. We also assume 
that the VEVs of all other bilinears are zero HgS], |2|. Therefore < Mf'''"^ >= Utafd^'f^. 
Following the above procedure for each term in the action (^) we have: 



1 X [^f^(-7.)^'^^-i,<,(7.)'''^-'1, (39) 



t I I 'I 1+^ 1 + ^ l+fj. t+fj, I fj, //i 1^ I fj, 5 



(40) 



< 0{AA) > 
< 0{AA)^ > 




0, 



< - 



2880 



29 



[tr{A,t.A^^)?tr[iA,t^A,^)^] > 
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< -^trliAM^ > = -^tUl, . (41) 

One can now change path-integration variables, from fermion to meson fields Ai. 
An important role in the dynamics of the system is played by the Jacobian of such a 




transformation, which is calculated in Appendix 2 following standard arguments ^6 
Including the Jacobian we get 

Zo = I[JdM.exp\-Y.i-^U^Ul^-^U^Ul^)-Y.2logU^ 

i y i,n i 

= UjdM.^^pi^-Y.i-lu^Ul^~^UtuU, + l\ogUtuU^)\, (42) 

which describes the dynamics in the (^2 = Pi = region of the phase diagram of figure 
|l]. To get a simple saddle point effective potential we complexify the condensate Ui and 
write, following ref. ||2^, UiUi+^ = U"^. This is justified, since we assume that, in the 
general case, the condensate Ui should be the same for all odd sites and all even sites 
separately. In our path integral the radial part of the contour is irrelevant (see Appendix 
2) and we can choose Ueven = Ue^"^ and Uodd = Ue~^'^. As we discuss in Appendix 2, 
the minimum value of the effective potential occurs for (f = 0. The zeroth order effective 
potential is defined Veff = —Seff/vol: 

143 

V,ff = 8lnU-AU'-—U^ (43) 

which is plotted in figure 0. 

We observe that it has a local maximum, but, as explained in p6|, this still implies 
stability of the broken SU{2) vacuum, due to the special properties of the Jacobian 
associated with the transformation from the \l/ variables to the meson variables Ai. 
This is reviewed briefiy in Appendix 2. 

From figure |^ it is evident that there exists SU (2) symmetry breaking even for the 
csae of /?2 = 0. This implies that the critical line in the phase diagram of figure |l] passes 
through the origin. This is also the situation argued to characterize the statistical model 
of 0, which may describe high-temperature superconductivity. 



3.2.2 Strong Coupling Expansion of 

We now look at the strong coupling expansion of the SU{2) field, up to and including 
order 0(/?|), 



p p p g 
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Figure 3: The effective potential for /?i = /?2 = 0. 



The zeroth order term has been calculated in the previous subsection (^). 
The first order term in the SU{2) integral is written 

= -f32N,l[JdV,,,I^'-{2^^) = -(32NpZo, (45) 

because the addition of a single plaquette will give an odd number of group elements on 
each side of the plaquette, and therefore will integrate to zero. 

At the second order the calculation is no longer so simple, 

^2 = f n / dV,,Xi^Vy^,) E(l - trV,) E(l - tr\/t), (46) 

by the same argument as used in the first order case, we can only have a non-zero integral 
where we avoid integrals over different numbers of Vs and V^^s. Thus the product Z]p(l — 
trVp) Eg(l - trVq) can be replaced by (N!^ - EptrVjtrVp) 

Z2 = ^ufdv.,Xi^Vy^,)iK + J:trKtrv;), 
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(47) 



Hence, the problem is to calculate the group integrals which make up the non-trivial 
part of Z2 namely; 



^2 = n / dv,,Xi'^Vy^^)T.tryptrv; 

i,fl p 
T.W j dV,X{2^) n / dV.,,lt\2^)trV,trVi 



1,1^ j 



(48) 



The full path integral is then written as 

We will be interested in the logarithm of this function: 



5, 



log 



zA\- /32iv„ + 



which can be expanded, ignoring constant factors and keeping terms up to 0(/?^ 



^e// = log ^0 + IJ- • 



(49) 



(50) 



(51) 



Therefore, the strong coupling expansion leads us to calculate the group integrals around 
the plaquette in the function Jp 



(52) 



n.,^ep/rf^i,M^^''(2V^) 

The denominator has already been given earlier and is the product around each side of 
the plaquette of the {AA) polynomial in (pBI). For convenience let us label the plaquette 
as having sides 1,2,3,4 and sites A,B,C,D 
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(1 + -trA^A^ + • ■ ^^'(l + -trA^A^ + • ■ x 

x(l + UrA^M + ■ ■ ■)"'(! + frA^Ai + ■ ■ ■)"\(53) 

These brackets can be expanded and truncated to 0{AAY since they still contain the 
grassmann fields (\E'\E'). The numerator contains a multitude of group integrals which 
need to be evaluated around the plaquette. The method, outlined briefiy in Appendix 3, 
is somewhat involved and the interested reader can refer to for the algorithms used. 
Suffices to say that the terms which have an even number of (\&\E') are kept since they 
will form the meson states. 

After a tedious calculation along the above sketched lines, the quantity Jp becomes 

= 1 + ^^{UaUbUcUd? + '^{UAUBUcUofiUlUl + U^U^ + U^Ul + UlUl) 

+ ^-^{UaUbUcUd)\ (54) 

and since in three dimensions we have one plaquette per link the effective potential to 
0(/5|) is (again with UiUi^^ = f/^): 

This is plotted in ffg. |^, with (32 taking a range of values between and 0.5. The 
behaviour does not change qualitatively as we increase (32-, showing that the symmetry 
remains broken as we move up the (32 of ffgure |[ This is as expected assuming a continuous 
critical line. 

3.3 The Phase Diagram for = 0, A 7^ 

Let us now complete our analysis on the phase diagram by concentrating on the region 
of strong SU{2), P2 = 0, keeping Us{l) coupling arbitrary (bottom horizontal axis of ffg. 
1). In this part of the phase diagram one can integrate out the (strongly coupled) SU{2) 
gauge ffelds to derive an effective action for the fermion and Us{l) gauge ffelds. The 
SU{2) path integration is performed along the lines of ref. |2^. In the strong coupling 
limit for SU{2), P2 = 0, the effective action, obtained after integration of the SU (2) gauge 
ffelds, reduces to the a sum of one-link contributions, S^fj = Seff{A, A), with 

A^{x)t = ^b{x + a)-f^Ul^^''ix) Mx)t = ^6(x)(-7^)t/,,^^"(x + a) (56) 
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Figure 4: The effective potential for /3i = and for (32 = 0,0.125,0.25,0.375,0.5, the 
corresponding curves lying in order of decreasing magnitude of their maxima from right 
to left in the figure. 



where U^,^ denotes the Us{^) group element, a is the lattice spacing, and the latin indices 
a, h are colour SU{2) indices. Below we shall proceed to evaluate explicitly this strongly- 
coupled effective action along the /32 = axis of the phase diagram in fig. |I|. 

For the SU{N) case the effective action exp{—Seff) is known in an expansion over A, 



A [26|. This will be sufficient for our purposes here: 



Seff = ^Tr{AA) + ^^^^^ _ ^^ [-Tr[{AAf] + ^{Tr[AA]f] + . . . + ^{detA+detA) + . . . 

(57) 

The determinant terms are associated with baryonic states p6|. We also note that for 



the U{N) case the determinant terms are absent. In the phase diagram of fig. 1 the 
U{2) case occurs at the point /?2 — ^ 0,/?i 0. We approach this point asymptotically, 
by working on the (32 = line, and assuming (3i ^ 0. We first notice that the Abelian 
phase factors of the ^/^(l) interactions cancel from the expressions for the traces of ^, A 
in the effective action (ISTf ). Moreover, from the discussion of section 2, we know that the 
SU{2) (strong-coupling) integration cannot produce a parity-invariant condensate, since 
the latter is not an SU{2) singlet p|. The resulting effective action should be expressible 
in terms of SU{2) invariant fields. Thus, on the axis (32 = there is no possibility for the 
Us{l) group to generate a fermion condensate. This implies that for very strong SU{2) 
group the symmetry is restored for arbitrary Us{l) couplings. 

This is a very important fact, indicating the existence of a not well-defined limit 
P2, Pi 0, since from the discussion in this and the previous subsections it seems that 
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there is an ordering problem in how one approaches the point [Pi, (32) = (0,0). This 
indicates that the shape of the critical line around that point is the one depicted in fig. 
|l], concaving upwards. By continuity arguments, then, one expects the shape of the entire 
critical line to be the one depicted in the figure. About the point {Pi, P2) = {Pi, 0), where 
the SU{2) interactions are negligibly weak, and thus irrelevant, one expects an almost 
vertical shape of the critical line. 

This discussion completes our analytical results for the phase diagram of the SU{2) 
Us{l) gauge theory. As we have already mentioned above, the derivation of the precise 
shape of the critical line, separating the phases of unbroken SU{2) symmetry from the 
region where symmetry breaking occurs, requires a proper lattice simulation analysis, by 
means of a fermionic algorithm. We hope to be able to address these issues in a future 
publication. However, the above results will be sufficient for our purposes in this work, 
and will enable us to present physically intreresting scenaria, pertinent to the physics of 
high-temperature superconductivity, which we shall discuss in section 6. 

Before doing so, it will be essential to review in the next two sections : (i) the symmetry 
breaking properties of our non-Abelian gauge model from the point of view of Goldstone's 
theorem and the existence of a local order parameter, and (ii) the role of non perturbative 
effects. These will be important in considering the coupling of the model to external 
electromagnetic fields, as required for the study of superconducting properties. 

4 Kosterlitz-Thouless Realization of Superconductiv- 
ity in the SU{2) Us{l) model 

This section is mainly a review of results that appear in the literature regarding the 
model P, ^ |lT). It mainly serves as a comprehensive account, for the benefit of the non- 
expert in the area, of the various delicate issues involved, which play a very crucial role 
in the underlying physics. 

An important issue in the model (§) is the existence of a global conserved symmetry, 
namely the fermion number, which is due to the electric charge of the fermions ^ . The 
corresponding current is given by 

J, = j2'^'l,^c (58) 

c=l 

This current generates a global Ue{^) symmetry, which after coupling with external elec- 
tromagnetic fields is gauged. 

In the absence of such external potentials, the symmetry Ue{^) is broken spontaneously 
in the massive phase for the fermions ^ . This can be readily seen by considering the 
following matrix element (see figure §): 

5" =< 5;|J,|0 >, a = 1,2,3 ; J^, = ^7/.^ (59) 



22 



Figure 5: Anomalous one-loop Feynman matrix element, leading to a Kosterlitz-Thouless- 
like breaking of the electromagnetic t/em(l) symmetry, and thus superconductivity, once 
a fermion mass gap opens up. The wavy line represents the SU{2) gauge boson B^, 
which remains massless, while the blob denotes an insertion of the fermion- number current 
= \1/7^\1/. Continuous lines represent fermions. 



As a result of the colour group structure only the massless gauge boson of the SU{2) 
group, corresponding to the (T3 generator in two-component notation, contributes to the 



graph. The result is |11 



S =< Bl\J,\0 >= (sgnM)e^,p^ (60) 



where M is the parity-conserving fermion mass (or the holon condensate in the context of 
the doped antiferromagnet). In our case this mass is generated dynamically by means of 
the Us{l) interactions, as we discussed above, provided the coupling constants were lying 
in the appropriate (strong) regime of the phase diagram of gig. |l[ 

The result ( |60D is exact in perturbation theory, in the sense that the only modifications 
coming from higher loops would be a multpilicative factor i_li(^p-^ on the right hand side, 
with n(p) the i?j^-gauge-boson vacuum polarisation function |]11 . 

As discussed in [|Hll|, the B^ colour component plays the role of the Goldstone boson oi 
the spontaneously broken fermion-number symmetry. If this symmetry is exact, then the 
gauge boson B^^ remains massless. This is crucial for the superconducting properties [Q, 
given that this leads to the appearance of a massless pole in the electric-current two-point 
correlators, the relevant graph being depicted in figure 

It can be shown [Q that in the massive-fermion (broken SU{2)) phase, the effective 
low-energy theory obtained after integrating out the massive fermionic degrees of freedom 
assumes the standard London action for superconductivity, the massless excitation being 
defined to be the dual of B^^: 

d^(p = e^,,pd,Bl (61) 
All the standard properties of superconductivity, Meissner effect (strongly type II ||^), 
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Figure 6: The lowest-order contribution to the electric current- current correlator 
{^\J ij,{p) J v{~P)\^) ■ The blob in the propagator for the gauge boson indicates fermion 
loop (resummed) corrections. The blob in each fermion loop indicates an insertion of the 
current J^. 

flux quantization and infinite conductivity, follow then in a standard way, provided the 
excitation (and, hence, B^) is exactly massless. 

Having discussed the spontaneous breaking of fermion number symmetry in the mas- 
sive fermion phase, it is natural to enquire on the nature of symmetry breaking, in the 
sense of establishing the existence or absence of a local order parameter. In this re- 
spect, our discussion will parallel that of ref. @]. The neutral parity-invariant condensate 
< — \l/2\l/2 >5 generated by the strong Us{l) interaction, is invariant under the 

U{1) ® Ue{1), as a result of the coupling of Bf^ in the action, and hence does not 
constitute an order parameter. This is a characteristic feature of our gauge interactions. 
Putative charge 2e or — 2e order parameters, like the pairing interactions among opposite 
spins in the statistical model of @, ^, e.g. < \l/i\l/2 >, < ^1^2 > Q will vanish at any 
finite temperature, in the sense that strong phase fluctuations will destroy the vacuum 
expectation values of the respective operators, due to the Mermin- Wagner theorem. Even 
at zero temperatures, however, such vevs yield zero result to any order in perturbation 
theory trivially, due to the fact that in the context of the effective B^ gauge theory of 
the broken SU{2) phase, the gauge interactions preserve 'flavour'. For a more detailed 
discussion on the symmetry breaking patterns of (2 -|- l)-dimensional gauge theories, and 



the proper definition of order parameter fields, we refer the reader to the literature |TT| , g . 

Thus, from the above analysis it becomes clear that gap formation, pairing and su- 
perconductivity can occur in the above model without implying any phase coherence. 



^In four-cmponent notation, such fermionic bilinears correspond to < ^'75 5" >, < ^754" >, considered 
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5 Instantons and the fate of Superconductivity 



An important feature of our model is that, due to the non-Abehan symmetry breaking 
pattern SU{2) — > f/(l), the abehan subgroup U{1) G SU{2), generated by the Pauh 
generator of SU{2), is compact, and may contain instantons [|T^, which in three space- 
time dimensions are hke monopoles, and are known to be responsible for giving a small 
but non-zero mass to the gauge boson B^, 

~ e"^^° (62) 
where 5*0 is the one-instanton action, in a dilute gas approximation. Its dependence on 



the coupling constant g2 = gsu{2) is well known |T9| : 



^0 (63) 

For weak coupling g2 the induced gauge-boson mass can be very small. However, even 
such a small mass is sufficient to destroy superconductivity, since in that case there is no 
massless pole in the electric current-current correlator. 

The presence of massless fermions, with zero modes around the instanton configu- 



ration, is known |]T9| to suppress the instanton effects on the mass of the photon, and 
under certain circumstances, to be specified below, the Abelian-gauge boson may remain 
exactly massless even in the presence of non-perturhative effects, thus leading to super- 
conductivity, in the context of our model. This may happen [jl9l if there are extra global 
symmetries in the theory, whose currents connect the vacuum to the one -gauge-boson 
state, and thus they break spontaneously. This is precisely the case of the fermion number 
symmetry considered above. In such a case, the massless gauge boson is the Goldstone 
boson of the (non-perturbatively) spontaneously broken symmetry. 

However, in our SU{2) ® Us{l) theory, discussed in this work, as a result of the 
(infinitely strong) Us{^) interaction, a mass for the fermions is generated, so we are not 
facing a problem with zero modes. Our analysis is based on a Wilsonian treatment, 
where massive degrees of freedom are integrated out in the path integral. This includes 
the gapful fermions, and tha massive SU{2) gauge bosons. The resulting theory, then, is 
a pure gauge theory U{1) G SU{2), and the instanton contributions to the mass of 
are present, given by (^), in the one-instanton case. 

We now remark that Supersymmetry is known to suppress instanton contributions. 
For instance, in certain = 1 supersymmetric models with massless fermions, considered 
in ref. the instanton-induced mass of the Abelian gauge boson is given by: 

^^gauge boson ~ 6 (64) 

which is suppressed, compared to the non-supersymmetric case (^). 
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N = 2 supersjTiimetric theories in three space-time dimensions constitute additional 
examples of theories where the abelian gauge boson remains exactly massless, in the 
presence of instantons ||19|, |39|. Such theories have complex representation for fermions, 
and hence are characterized by extra global symmetries (like fermion number). In view of 
our discussion above, such models will then lead to Kosterlitz-Thouless superconductivity 
upon gauging the fermion number symmetry. 

We also remark that in supersymmetric theories of the type considered here and in 
ref. it is known that supersymmetry cannot be broken, due to the fact that 

the Witten index (—1)^, where F is the fermion number, is always non zero. Thus, 
in supersymmetric theories the presence of instantons should give a small mass, if at 
all, in both the gauge boson and the associated gaugino. However, in three dimensional 
supersymmetric gauge theories it is possible that supesymmetry is broken by having the 
system in a 'false' vacuum, where the gauge boson remains massless, even in the presence 
of non perturbative configurations, while the gaugino acquires a small mass, through non 
perturbative effects. The life time, however, of this false vacuum is very long and 
hence superconductivity can occur, in the sense that the system will remain in that false 
vacuum for a very long period of time, longer than any other time scale in the problem. 

The short reviews of symmetry-breaking patterns and the role of non-perturbative ef- 
fects, just presented, provided us with the necessary equipment to attempt a construction 
of possible scenaria, which might simulate the interesting physics underlying the high- 
temperature superconducting cuprates. A rather preliminary and heuristic discussion 
will be presented in the next section. A more detailed analysis, especially in the context 
of the statistical models of ref. , requires proper lattice simulations which automatically 
incorporate non perturbative effects. This, however, falls beyond the scope of the present 
paper. 



6 Application to the Physics of high-temperature su- 
perconductors 

6.1 Phenomenology of high-temperature superconducting ma- 
terials 

In this section we would like to consider a possible application of the above SU{2) 
Us{l) model H] to the physics of high-temperature superconducting cuprates. Recent 



experiments ||2^ have demonstrated an extremely unconventional and rich structure of 
these materials, not in their superconducting phases, but rather in the normal phase. The 
phenomenology of the high-temperature cuprates may be summarized by the temperature- 
doping concentration phase diagram, shown in figure |^. 

The phase diagram shows clearly a very-low (including zero) doping antiferromagnetic 
phase (AF). Above a critical doping concetration (point A in fig. ^, AF order is destroyed. 
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Figure 7: The temperature- doping phase diagram summarizes the current (experimentally 
observed) situation in high-temperature superconducting cuprates. Notice the existence of 
an intermediate zero-temperature phase, characterised by the existence of preformed pairs, 
leading to a pseudogap. 



but the interesting issue is the existence of a phase, named 'psuedo-gap phase', which 
interpolates between the AF and the sperconducting phases (dSC), the latter being known 
to be of (i-wave type [|1^ . 



It is a general belief today, supported by many experimental results [^, e.g. results 
on optical conductivity, photo-emission, transport etc., that the pseudo-gap phase is char- 
acterised by pairing ('pre-formed pairs'), leading to the existence of a mass (pseudo)-gap 
in the fermionic spectrum, which however is not accompanied by phase coherence. This 
situation is in sharp contradction with in standard BCS theories of superconductivity, 
according to which phase coherence appears simultaneously with the appearance of a 
gap. 

The pseudogap phase is separated by a critical temperature curve Tc (doping) from the 
d-wave superconducting state, characterized by the sharp drop in resistivity, and strong- 
type II superconductivity (penetration depth of external magnetic fields is of order of 
a few thousands of Angstroms). Today, the general belief is that the superconducting 
pairing is of BCS type involving four-fermion itneractions among the charged excitations. 
However, the four-fermion interactions do not have to be phononic. 

The pseudo-gap phase is also separated by another curve T* (doping) from the normal 
state phase, where there is no gap, but where there are abnormal normal state properties, 
such as linear dependence of the electrical resistivity with temperature for a wide range 
of temeprature scales etc. All such properties point towards a non-fermi liquid behaviour 
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of the normal state, which is experimentally observed, as far as we understand, at least 
in the regime of optimum doping (shown in the figure |^). 



6.2 Strongly-coupled U (2) gauge theory and the pseudogap phase 

In this section we shall argue that the gauge theory SU{2) (g) Us{l) of ref. 0, whose 
low-energy limit has been studied in this paper in some detail, may provide a satisfactory 
qualitative explanation of the phase diagram of fig. |^, especially as far as the appearance 
of a pseudogap phase is concerned. For the purposes of this article, we shall concentrate 
in the zero temperature region of the graph. Our method will be that of ref. i.e. 
approaching the pseudogap phase by studying the excitations about the nodes of the d- 
wave superconducting gap. We shall not deal here with excitations away from the nodes 
of the gap. Our hope will be that similar (long range) gauge-interaction phenomena are 
responsible for the formation of the bulk of the d-wave gap and the pseudogap. From the 
preliminary finite-temperature analysis of ref. it becomes clear that the gaps that 



open up at the nodes disappear at much lower temperatures (T < O.IK) than the bulk of 
the (i-wave type gap (Tc = O{100K)), and this means that the predictions made in this 
work, if true, can be realised only if one looks at low temperatures 

Before starting our analysis on the (zero-temperature) pseudogap phase, we should 
point out that the gauge theory at hand, is also in agreement with deviation from fermi- 
liquid behaviour in the normal (no mass gap) phase. Indeed, as we discussed in previous 



articles a U{1) fermion-gauge theory in (2 + 1) dimensions, in which the mass of the 
fermions is generated only dynamically through the gauge interactions, is characterised 
by non-trivial infrared fixed points, which according to general arguments |Q is sufficient 
to drive the theory away from the Landau fermi liquid (trivial infrared structure). 

Let us continue our discussion on the diagram of fig. by considering the Antiferro- 
magnetic phase (AF). In such a phase, the only excitations are assumed to be spin degrees 
of freedom. There are no charged excitations, and the pertinent dynamics is described 
by the magnon z sector of the model (^. As the doping exceeds a critical concentration 
(point A in fig. |^), the antiferromagnetic order is destroyed. In the context of simple low- 
energy CP^ models, which describe adequately the dynamics of the AF sector, this can be 
seen easily by applying renormalization-group arguments ^ , and taking into account 
the dependence of the respective coupling constant on doping, in the way explained in 
refs. i, 0. 

The important question is whether superconductivity does not set in immediately, but 
one has to pass through the intermediate phase AB, where a 'pseudo-gap' appears, but 
no phase coherence exists. As we shall argue now, our strongly-coupled gauge theory 
SU{2) Us{l), presented above, may offer an explanation for the phenomenon. 

^However, we point out that the presence of external magnetic fields may enhance these values pof , 
as, for instance, is the case of the experiments involving thermal conductivity measurements See 
discussion at the end of this section. 
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To this end, we first remark that above the critical doping concentration that marks the 
on-set of disorder [A in fig. the z magnons are massive, with masses which themselves 
depend on the doping concentration |^5|, |15|, In this regime, there are both charge and 
spin [z field) excitations. Integrating out the massive magnons z, the long-wavelength 
dynamics of the charge excitations is described by the effective SU (2)®?7s(l) gauge theory 
of ref. (§P- The gauge Us{l) interaction is capable of inducing dynamical opening of 
a holon gap (pairing) if the pertinent coupling constant of the statistical model lies inside 
the SU (2) broken regime of the phase diagram of fig. |T[ 

We now remark that in the statistical model (|^), which will be the basis of our dis- 
cussion in this section, the inverse couplings of the SU{2) and Us{l) gauge groups he in 
the straight line AB depicted in fig. ||, as a result of (^. In the condensed- matter model 
of ref. 0, then, the local gauge group is U{2) rather than SU{2) ® Us{l) Q. 




Figure 8: Phase diagram for the SU{2) x Us{l) gauge theory, viewed as a low-energy 
continuum limit of the solid-state model for doped antiferromagnets of ref. [6J. The 
straight line indicates the specific relation of the coupling constants in the model. 

At present, the precise shape of the critical line is not known, since it requires the 
construction of an appropriate fermionic algorithm, which will allow for a proper lattice 
study of the model. The strong coupling analysis in this paper has demonstrated, however, 
that the critical line passes through the origin of the graph, concaving upwards in the 
way shown in the figure. Also, from the behaviour of the critical line about the point 
{P2 = 00, Pi), it is evident from the graph of figure ^ that the intersection point B defines 

^This, however, does not affect the resuhs of the previous analysis pertaining to the mass generation. 
The only difference of the U{2) case is the absence of haryons from the spectrum [ p6| . 
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an upper bound for the inverse coupling in order for the system to be in the mass- 
generation phase. For all practical purposes it is qualtitative meaningful to assume an 
almost vertical shape of the critical line at the intersection point B, which implies that 
the critical coupling for mass generation for the coupling jSi in the statistical model is 
still given by the single U{1) gauge theory critical coupling, i.e. /?i < jS^ (see figure |I]). It 
is known that ~ 32/7r^ so, on account of our discussion in subsection 3.1 and 

d^), such a gauge pairing would occur in the following range of doping concnetrations: 

6af<6< = 1-^7 (65) 

where 5af denotes the doping concentration at which the AF order is destroyed. 

This phase is characterised by the breaking of chiral symmetry. However, as discussed 
in IjTT], Q , and reviewed in section 4, the symmetry breaking occurs without a local order 
parameter. Strong phase fluctuations destroy the putative order parameter for 2 + 1- 
dimensional QED. This is an exclusive feature of the 2 + 1 gauge interactions, and as 
we argue now, it is responsible for the appearance of a pseudo gap. The concept of 
the pseudogap is associated precisely with the presence of a non-vanishing mas gap and 
thus the existence of 'pre-formed pairs' in underdoped cuprates ^Mj] , but in absence of a 
local order parameter (phase coherence) in the model. The situation is analogous to the 



Kosterlitz-Thouless mode of symmetry breaking |12|. The important issue to understand 



is why there is no superconductivity in the model. 

This issue is related to the presence of instantons, discussed in the previous section, 
which are responsible for giving a small but finite mass to the gauge boson B^, and thus 
destroying the basic criterion of superconductivity (see fig. |^). The presence of fermions 
does not change this. In the broken SU{2) (gapped) phase of the model (^ the fermions 
are already massive, due to the extra Us{l) interactions, and hence there is no issue of 
zero modes that could screen the instanton effects. From (|62D, as well as the fact that the 
one-instanton action exhibits the following dependance on the gsu(2) coupling constant 
for the case at hand |T^: 

^o~^^^^ ; ^1^(2) oc J(l -5) (66) 
9su{2) 

one observes that the instanton-induced -Bj^-boson mass decreases upon increasing the 
doping concentration 5 in the sample. 

In the context of the statistical models of ref. 0, etc, one should also consider the 
coupling of superconducting planes, by means of Coulomb interactions among the charge 
carriers (electrons). Such interactions may result in a small leakage of electrons across the 
planes, which inevitably leads to fermion-number non-conservation on the plane. In ref. 
|28| , within a spin-charge separating framework, such an interplanar coupling has been 
represented by inserting in the path-integral a term of the form: 

j t^^e2'/'^'^^*"^-+'^-=- (67) 
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where a = 1,2 runs over 'colours' in the model of SU{2), and 1] is a. Majorana spinor. Due 
to this, the term ( |67D in the effective action violates fermion number, and is interpreted 
as implying a hopping of both spin (z) and charge (\E') degrees of freedom. Notably, t] 
may play the role of the supersymmetric partner of the gauge boson, in a = 1 
supersymmetric formulation of the model 0, which is possible upon certain relation ||28[ 
among the couplings of the microscopic spin-charge separating model of 

The explicit breaking of the fermion number symmetry by the interplanar coupling, as 
well as the absence of fermion zero modes in the massive phase (due to the Us{l) interac- 
tions) imply that the presence of fermions will not cancel the instanton-induced small mass 
of the gauge boson (^) [|. In such a case, then, the gapped phase will be characterised by 
the presence of pairing, mass gap, but no phase coherence and superconductivity, features 
shared by the pseudogap phase observed in cuprates (fig. |^). 

The above considerations are rather heuristic at present. The complete analysis would 
necessitate a lattice simulation of the model (0) in the presence of instanton effects (in 
the broken phase after mass generation due to Us{l) interactions). Analytical results at 



present exist only for A^ = 2 supersymmetric theories, as we mentioned above [O, 39 



which however seem not to correspond to the physics of the cuprates which appear 
to have at most A^ = 1 supersymmetry upon coupling of the superconducting planes. 
Supersymmetry supresses instanton effects in some cases (0), and some times may lead 
to a massless gauge boson, although such a case at present seems to characterise N = 2 
supersymmetric theories. We hope to return to a more detailed study of such issues in 
the future. 

The suppression of instanton effects by supersymmetry, which in our class of statistical 
models may occur for certain doping concentrations points to the following possibility: 
As one increases the doping concentration, a region is reached where there is a special 
relation among the various coupling constant of the effective spin-charge separating theory. 



leading to a A^ = 1 supersymmetry p8[ For instance, in the context of models of ref. [Tl 
such a supersymmetric point could be reached for doping concentrations 6*, such that 
t' ~ v07'(l - 5)3/^ where the prime denotes next-to-nerarest neighbor hopping (t) and 
Heiseneberg exchange energies (J). By tuning the couplings one may arrange- always in 
the cotnext of phenomenological models- for a situation in which 6* < 5^^^\ This would 
imply that, within the region of dopings for which the gauge statistical interactions are 
responsible for the opening of a gap and pairing, the suppression of instanton effects due to 
the presence of (supersymmetric) fermions may be sufficient to allow for a gauge-theory- 
induced Kosterlitz-Thouless (KT) superconducting gap at the d-wave nodes. The KT 
nature of the gap implies that once opened such a gap cannot affect the ci-wave character. 



^This supersymmetry carries non-trivial dynamical information about the spin-charge separation 
mechanism underlying the model, and hence it is different from the non-dynamical global supersym- 
metry algebras, at specific points of the coupling constants, discovered in |Q. 

^Although, a reduction of order (^) might be expected in = 1 supersymmetric cases [p8| , occuring 
for particular values of doping. See discussion below. 
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This scenario for superconductivity has been advocated in ref. 0. In a related, but less 
probable, scenario, the tunnelling to a 'false' supersymmetry-broken vacuum |T^ could 
occur in the dSC region of the phase diagram (see fig. 0), 



6.3 Superconducting phase and additional four-fermion interac- 
tions 

Despite these appealing scenaria for the role of gauge (spin-spin) interactions for inducing 
KT superconducting gaps at the nodes, in the realistic situation the onset of ((i-wave) 
superconductivity occurs at higher doping concentrations, for which the attractive four- 
fermi couplings among charged excitations in the effective fiedl theory become strong 
enough, so as to overcome the gauge interactions, and lead to a standard BCS type 
pairing among the charged excitations. 

In such a case one then is forced to consider the effect of additional contact interactions, 
among the holons, which are up and above the gauge interactions considered so far. 
The simplest, and most likely the most relevant, of such interactions are contact four 
fermion interactions ^. From a microscopic point of view such interactions may 
arise as an effective way to describe the tendency of holes to occupy nearest neighbor 
lattice sites ^ |l8l, or in some recent scenaria they may describe attractions due to 
screened Coulombic interactions among charged excitations, as a result of the interplanar 



coupling [^. The corresponding coupling constant will again depend on the doping 
concentration and this will imply interesting phase structures. 

The additional four-fermion interactions are, then, viewed as being responsible for the 
appearance of an ordered c?-wave state, with the gap being characterized by four nodes. 
We shall be fairly phenomenological, given that a detailed incorporation of extra four- 
fermion interactions in our strongly-coupled gauge model is not striaghtforward. However, 
a phenomenological analysis will be sufficient to demonstrate the main unconventional 
features of our scenario for high-temperature superconductivity, at least in the context of 
a continuum effective field theory. 

It is known that, in the context of relativistic models we consider here, as a result of 
linearization about the nodes of a d-wave superconducting gap, the four-fermion interac- 
tions are the only ones which become renormalizable (relevant) in the framework, 
where is a flavor number for fermions. As an instructive example, consider, for instance. 



Gross-Neveu type four-fermion couplings in the effective lagrangian pi 

o 

/ \ 2 



^4/ = «:E(^a*i (68) 



a=l 



where \E'a are the relativistic spinors (^ describing the excitations about the nodes of a 
d-waye gap. 

From standard arguments on the phase structure of Gross-neveu type couplings. 



we are considering here, we know that pair formation, and hence mass generation, in 
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four-fermion theories occurs for dimensionless inverse couplings, A = weaker than a 
critical value, 2/7r^. However, the full phase diagram, incorporating the SU{2) and Us{l) 
couplings as well, as appopriate for the model of will be more complicated. However, 
for our purposes in the present work it will be sufficient to consider only the effects of the 
Us{l) coupling, responsible for the mass generation in the model. 

A phase diagram for an abelian gauge theory with extra four-fermion interactions (of 
Gorss-Neveu type) has been derived analytically in the context of a Schwinger- Dyson large 
fermion- flavour analysis in ref. [Q, and in principle the result can be checked in lattice 
models. For our purposes in this work it will be sufficient to assume the validity of the 
large- flavour- number continuous results of ref. 0] , and concentrate on the pertinent phase 
diagram in the coupling constant space between the gauge g and four-fermion couplings 
K, depicted in fig. ||. 




Figure 9: A generic phase diagram of the theory with Us{l) gauge and four-fermion 
(Gross-Neveu) interactions. The critical line separates the phase of unbroken symmetry 
from that of broken symmetry. The symmetry breaking is due to the fermion condensate. 
The dotted line is conjectural at present, and indicates the on-set of a local order parameter 
due to the dominance of four-fermion interactions. 



In toy models of doped antiferromagnets |I8|, which are sufficient for our illusttra- 



five purposes, such Gross-neveu four fermion terms are expressing the tendency of holes to 
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break as less bonds as possible in the antiferromagnetic lattice, which is the configuration 
featuring the holes sitting next to each other. Such interactions, may be described by 
adding to the Hamiltonian attractive four-fermion interactions of the form: 

-k:4Mj) ■.■.4Mj + '^) ■■ (69) 

where ipa,o: = 1,2 are Grassmann (holon) operators, j denote lattice site, and : ... : 
denotes normal ordering of quantum operators. The normal ordering conventions are 
such that a fermion bilinear is written as: 

4^ =: 4^ '■ + < i^H >, < >= ^ (70) 
with S the doping concentration in the sample. Such terms may be asssembled, in the 



continuum, low-energy, limit JTS], ^, into Gross- Neveu four-fermion terms of the form 
(|68|) , where the spinors are constructed as in (^. 

At this stage, the coupling constant k is a phenomenological parameter. However, 
from quite generic arguments, one would expect it to increase upon increasing the doping 
concentration 6 in the sample, since the larger the doping, the bigger the probability of 
the holes to lie in adjacent sites of the lattice. 

At present, the only case where four-ho Ion-operators appear with well-defined coupling 
constants in terms of the microscopic parameters of the theory, is the t — j or Hubbard 



model case, where, however, the four fermion interactions are repulsive [|T^, In the 



models of ref. |T^, for instance, such Hubbard four-fermion couplings HHubbard assume 
the generic form: 

^Hubbard OC Ko(t', J')- r (71) 

1 — 

where /to(i', J') is an appropriate function of the next-to-nearest hopping element, and 
Heisenberg exchange energies. 

Combining (^) with such repulsive interactions, one then arrives at a generic coupling 
for four-holon (Gross-Neveu) operators in the model of the form: 

^,f = f,{6)-f^o{t',J')-^ (72) 

1 — 

where k{6) is, at present, a phenomenological parameter, which however is expected to 
increase, as we said, with increasing S. For the four-fermion interactions to be attractive 
one needs K4/ > 0, and this places restrictions on the regime of doping, for which the 
interactions are going to lead to dynamical mass gap generation. 

When combined with the phase diagram of fig. ^ this implies that pairing due to 
four-fermion interactions would occur for doping concentrations in a region determined 
by the critical line of fig. 

By appropriately choosing in the context of phenomenological models, it is then 
possible to arrange for a situation like the one depicted in fig. |^, where the zero- 
temperature pseudogap phase interpolates between the AF and the standard BCS type 
(i-wave superconducting theory. 
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Notice that the dynamical mass generation due to four-fermi couphngs leads to a 
second order transition, at zero temperatures, and hence to phase coherence, as is standard 
in BCS type of pairing. This should be contrasted with the gauge situation described 
above, which leads to Kosterlitz-Thouless type of breaking [^, ^, It would probably 



imply the existence of a cross-over line (dotted) in the diagram of figure separating 
the region of the broken symmetry phase where a local order parameter is present, due 
to the dominance of the four-fermion itneractions, from the region where the Kosterlitz- 
Thouless mode of symmetry breaking (absence of a local order parameter) occurs, due 
to the dominance of the gauge interactions. Such phase diagrams should be confirmed 
by detailed lattice simulations, using appropriate fermionic algorithms, which fall beyond 
the scope of the present work. 

The interesting feature is that experimentally one can make a distinction between a 
gap induced by the gauge interactions, or by four fermi interactions, as a result of different 
scaling of the mass gap with an externally applied magnetic field. A suggestive experiment 
along these lines is that of ref. [^], measuring the behaviour of the thermal conductivity, 
in both the superconducting and 'pseudo-gap' phases. Details on such issues are discussed 
inrefs. [H, |Ol|l|, |2l|. 



7 Conclusions 

In this work we have described a strong coupling expansion for an SU{2) ® Us{l) gauge 
theory, in three dimensional space time. From the physical point of view, such models may 
serve either as a prototype for physical applications to the physics of excitations about 
nodes in a d-wave high-temperature superconductor , or - when formulated in Euclidean 
lattices - as describing high-temperature phases of four- dimensional gauge theories of the 
Early Universe. 

Our analysis has indicated a phase diagram of the form depicted in figure |l|. The 
analytical results obtained in the present article pertained to the strongly coupled SU (2) 
sector. We have shown the absence of a finite critical coupling for SU (2) symmetry break- 
ing. The breaking of SU{2) is induced by dynamical parity-conserving mass generation 
due to the Us{l) strongly-coupled sector. The shape of the critical hne, spearating the 
phases of broken SU{2) symmetry, is still conjectural, since it requires proper lattice 
simulations with dynamical fermions, which is under consideration at present. 

An important ingredient in our analysis, which was motivated from our condesed- 
matter ancestor models 0, was the use of "naive" Dirac spinors on the lattice, and not 
Wilson fermions. The latter are known to violate explicitly parity-symmetry breaking. 



due to the Wilson term. This may lead to spontaneous violation of parity symmetry 



31 1, and therefore to a completely different phase diagram, although the issue is still 



unsettled 



similar effect may be induced by external electromagnetic interactions, which from a condensed- 
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An interesting application of our strongly-coupled gauge theory SU{2) (g) Us{l) was 
argued to be provided by a possible explanation of the (zero-temperature) pseudo-gap 
phase between the antiferromagnetic and ci-wave superconducting phases of the high- 
temperature cuprates (see fig. |^). Due to the special symmetry breaking patterns, 
SU{2) — i> [/(I) in the phase where a fermion mass is generated dynamically by the Us{l) 
statistical interactions in the model 0, and the existence of instanton configurations in 
the compact U{1) G SU{2), a small mass for the U{1) gauge boson can be generated. 
Such a small mass, although does not prevent pairing, however it spoils superconductiv- 
ity, since it leads to the disappearance of the massless pole in the electric current two-point 
correlators. As explained in the text, in such a sector, the gauge theory is responsible only 
for the opening of a mass gap in the fermion spectrum. However, the mode of symmetry 
breaking is of Kosterlitz-Thouless type, and hence no order parameter exists, since strong 
phase fluctuations destroy it (the parity-invariant fermion mass gap in the model of refs. 
is not an order parameter, since it is invariant under the respective symmetries, and 
hence there is no contradiction with the Mermin- Wagner theorem). Such a KT mass gap 
may be viewed as a pseudo-gap. In our scenario above, such a gap owes its presence to 
relativistic fermions at specific points of the fermi surface (e.g. nodes of the d-wave gap), 
argued to play a crucial role in the pertinent phase. We stress once more, that the impor- 
tant point in our approach was the Kosterlitz-Thouless nature of the gauge symmetry 
breaking induced by gauge interactions in (2 -|- 1) dimensions |Tl|, which discriminates 



our approach from others |22, 23 



Our belief that spin-gauge interactions may a play a crucial role in the underdoped and 
normal phase of the high-Tc cuprates is strengthened by the abnormal properties of these 
phases including the explicit observation of phase separation in the so-called stripe 
phase, occuring for a particular doping concentration |^3[. In the ci-wave superconducting 



phase, four-fermion BCS-like pairing may indeed occur, although the attractive four- 



fermion interactions, most probably, are not due to phonons, but of electronic origin |^ 

In the presence of both types of interactions, gauge and four-fermion, the effective the- 
ory model presents interesting phases. One way to determine the origin of the dynamically- 
induced mass gap in the various pahses is to study the behaviour of the system under the 
influcence of external fields, like in the experiment of PT|. It is known |ED|, BDI, 121]] that the 



matter point of view are natural to consider. Recently, the effects of constant magnetic fields on the 
opening of a mass gap at the nodes of d-wave high-temperature superconductors have been considered 
in the experiments of Mm. Claims that this may induce, for strong enough fields, a change of state 



of the condesate into a parity- violating one, have been made |52 . Indeed, in the case of a constant 
external magnetic field, perpendicular to the spatial plane, one has an external source term violating 
parity and time-reversal symmetry. For strong enough source fields it is possible that a parity-violating 
condensate is magnetically induced. Such a phenomenon is at present a conjecture, which needs to 
be demonstrated analytically (via Schwinger-Dyson analysis) or on the lattice. We postpone such an 
issue for future investigations. We should mention, however, that recent preliminary lattice |Q or 
continuum analyses, in the presence of an external field, showed that the magnetic field enhances the 
parity-conserving condensate. 
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gauge- field induced mass gap scales differently with an applied magnetic field as compared 
to the gap induced by four-fermi Gross-Neveu type interactions. Such a scaling may be 
determined by studying the thermal conductivity in the presence of an externally applied 
magnetic field, as in the experiments of [^]. Details of this analysis appear in | |21| |. 

We are, of course, aware that the simple, effective gauge field theory analysis we 
have just presented, may not be sufficient to explain quantitatively the rich phenomenol- 
ogy of the high-temperature cuprates. We believe however, that it constitutes a step in 
this direction. Our hope is that, due to the simplicity and universality that underlies 
the superconducting models based on the gauge symemtry approach, our results cap- 
ture essential features of the physical mechanism(s) underlying various phases of high- 
temperature superconductors, in much the same way as the single phonon BCS theory 
describes adequately the complicated physics of phononic superconductors. Moreover, as 
particle theorists, we also find this exercise very interesting, since it may imply that cer- 
tain phenomena, characterizing the Physics of the Early Universe, may have interesting 
counterparts in condensed-matter physics, and in particular the high-temperature (anti- 
ferromagnetic) superconductors. In this respect, specific mention should be made again 
to the work of Volovik [^, who pursues the analogies between Particle and Condensed 
matter Physics, by suggesting solid-state experiments, involving superfiuid Helium, as 
possible laboratory experiments which might shed light to the Physics of the early stages 
of our Universe. In the same spirit, the rich, and unconventional, structure of the high- 
temperature cuprates, depicted in fig. ^ may also find interesting, and possibly new, 
applications to Particle Physics. 
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Appendix 1: Rules for Multiple 'Meson' Pairing 



When we come to calculate the meson pairs on a particular site there is a subtle com- 
plication if there are two pairs on the site, ie. four ^s and four ^s. Let us take as an 
example the function [tr(Ai^^^j_^)]^ which lies on the link + /i). Expanded out this 
gives 

[^f'^(-7M)'^^^?;',^:;;(7M)'^''*f"'] X [^•"(-7M)'^^^f^,^fi(7M)'^'^^-'"] (73) 

Before, when we just had a term which produced one pair of mesons on each site there 
was only one unique way of combining the (^^) at that site to make the meson. However 
now, as we see, there are three ways at each site. Writing down the fields at site i; 

^^'^^^'"^'''^^'^'^^/'''^•^ ''^i^j'^'^^^'" (74) 
these can be combined into pairs in three ways 

% % % % % % % % 

or X *f'''^f'*f''^f''' 

or ^f^f'^^f'''^^'" X ^f'^^f'^f^^^'^''^ (75) 

when faced with a choice of three possible meson states, we must assume that each is 
likely to happen with equal probability (in a sense they can be seen as quantum states 
with equal energy) and so we write the final combination state at site i as 

3^? ' 3 * I * « 2 I I I 

+ i*,''^^f'"*f'''^''"*^^^f'*f''^^", (76) 

or 

\^tr,[Mf^Mr]tr.[MrMr] + \tr,[Mf'^Mr]tr.[MrMr] + \tr,[Mf^Mr]tr.[MrMr]] 

(77) 

This must be contracted with the gamma matrices on the link and then with the meson 
states on site i + i-i, calculated in the same way. 

This procedure must also be followed when the four (^^) pairs on the site come from 
two functions on adjacent sides meeting at that site. 
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Appendix 2: The 'Meson' Jacobian 



In this appendix we shall calculate the Jacobian of the transformation from the fermion 
variables '^a,a, ^a,a, to the meson variables A^^^ = ^^^o- This change of variables implies 
the following transformation in the path integral: 

n / d^id^^ ^ n / dMi, (78) 

i i 

we adapt the method outlined in ||2^. Let our initial path integral be written; 

n / d^^d^ie'^^ff^^''^ = He^-^^^I^l J d'^id'^,e'''^^%=o, (79) 

i i 

the above fermionic integral can be evaluated; 

J d^d^e*"-^^^ = J rf^d^e*^* = detJ ; ^J^ = j'^ba/s^bp ^ ^gg^ 

where detJ means the product of the eigen-values of J, which is a 2 x 2 matrix in two 
spaces. However there is a problem defining exactly what we mean by the measure dM, 
since M lives in two spaces. However if we instead regard M as a 4 x 4 matrix we can use 
some simple results from group integration [53, ^ to define a measure. The following 
identity is true for a ?7(4) group integral 

dU-^—e''^'^^ = detj', (81) 
detRU 

where R is any positive-definite hermitian matrix. The matrix j' is chosen to be the 
equivalent matrix in 4 x 4 space to the J defined in (^) which lived in two 2x2 spaces. 
But the determinant of both are the same, detJ = detJ. It is then a simple matter of 
associating, in (p]), dU t— > dM and RU i— *• M and we get 

d'^d^e'''^^ = I dM^—e''-^^, (82) 
J detM ^ ' 

and our path integral becomes 

n / rfMie-S«^"'^^**^+^=^/ . (83) 

Our path integral, dMi^ over the 4x4 space which we have represented as a group 
integral, can be viewed as a multi-component generalisation of a contour integral. 

If Mj had been a complex number, which is the case of a gauge theory with 



Kogut-Susskind lattice fermions, we would have had p6 



. , . ^ . (84) 
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This contour integral can be evaluated by parameterizing z = i?e*^, which is of course a 
representation of the U{1) group. In such a case, this property of the Jacobian, expressible 
in terms of contour Cauchy integrals, can be used to infer stability of the broken vacuum 
in strongly-coupled gauge theories, despite the appearance of a local maximum in the 
potential. A detailed discussion may be found in ||26|, where we refer the reader for more 
details. Below we shall only concentrate to a description of the basic results, pertinent 
for our analysis in this work. 




Figure 10: Explanation on stability despite the existence of a local maximum in the effective 
potential. 

In our problem, M is a 4 x 4 matrix and thus has 16 degrees of freedom. These degrees 
of freedom are illustrated in the expansion of M ([T^) in terms of the "lengths" (the ^s and 



JFs) along the 16 "axes". However we were interested in symmetry breaking along just 
one of these "axes", o"3.1, with a "length" given by the complex number ^3. So we can 
view the important part of the integral dM as being in the complex plane where the only 
degree of freedom is ^3. In this respect we have complexified the V.E.V. of ^3 in order 
to apply the contour integration. Eventually, the effective potential will be minimized for 
real < ^3 >. We have thus a standard contour integral as above (^i]). 

To make things clearer it is useful [^] to add an explicit chiral-symmetry-breaking 
fermion mass term, ^^m'^ia^^i, to our lagrangian. The VEV of Az{i) can be written 
as 



< A^ii) > = f/e*'^, even sites, 

< ^3(i) > = Ue-'"^, odd sites, (85) 

The VEV of the mass term < m\E'jCr3\E'j > is m < A3{i) > and so summed over even and 
odd sites will pick out the real part of the complex number: 
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'^m < A^ii) >= -mU COS ip . (86) 
i i 2 

Our effective potential now with the mass term is given in the form 

2 mU 
V^ff ~ -InU' - P{U') - cosy., (87) 

where P{U^) represents the polynomial in which makes up the rest of K//? given in 
the text (^3D,(|55|). It is now clear how the apparent maximum at the stationary point is 
interpreted. Since the path integral is effectively a contour integral, we can choose our 
contour around the circle \ < A3{i) > \ = U . 

Then, the important parameter for minimizing Ve/f is (f, and the minimum along the 
contour occurs at = 0. The local maximum in Ve/f lies along the radial direction U (see 
figure |lOD, which is irrelevent given that our variation of the potential is contrained to 



lie along the contour. As m ^ 0, the minimum flattens and the whole contour becomes 
degenerate. By this argument, one will have a dynamically generated non-zero VEV for 
^3, even in the absence of a bare mass term. Although this could equally be applied 
to any of the terms in the bi-linear expansion of Mi, we know from the discussion of 
Vafa-Witten [^, briefly reviewed in section 2, that the parity-conserving mass ^3 is 
energetically favoured. 
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Appendix 3: Outline of Strong-Coupling Computa- 
tional Rules 



We want to calculate the following function on the plaquette 



Jp 



with 



2 3 4 



4 36 576 

Our convention will be to label the sides of the plaquette 1,2,3,4 and the sites A,B,C,D. 
We need to perform the group integrals / dVi^^ on each side of the plaquette. In the case 
of this will give a product of polynomials in 0{AA) on each side 



Zo= ni,^ep/t^KX(2v^) = (l + ^tr(AiAi) + ---)(l + ^tr(AiAi) + ---) X 

x(l + Ur{A^A:,) + ■ ■ ■)(! + \tr{A^A^) + ■ ■ ■) (90) 

Zp will not be separable into products of simple link-polynomials becasue of the trVptrVj 
which connects the colour indices of the links around the plaquette in a non-trivial way. 
For example there could be a term like tr{AiA2A3A4A4^A3A2Ai) in Zp which would defi- 
nitely be absent in Z^. 

Zp and Zp are both truncated due to the Grassmann fields (\1/\E') contained in the 
A, A functions. So eg. a term like tr{AiAi)Hr{A2A2) will be immediately zero because 
there are too many fields at site B between sides 1 and 2. However eg. a term like 
tr{AiAi)Hr{AsA3)^ will survive. 

Alot of the integrals in Zp will turn out to be equal to their equivalent values in Zp. 
In fact the group integrals without at least a term of first order in (AA) on every side will 
just be equal to their zeroth order result. This is obvious if we consider eg. 



dVidV2dV3dV4{tr{AV^)tr{AV))hrVptrV^ or o(aa) 



(91) 



The diagram represents the plaquette with the clockwise flowing arrows being trVp and 
the anti-clockwise trV^ , (the reader should refer to Creutz for a description of dia- 
grammatic group integration). We do the group integral on one of the blank sides first. 
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Following Creutz [54|, this produces delta functions which contract with the other group 
matrices in trVptrV^ all the way around the plaquette, removing them. This leaves the 
initial {tr{AV^)tr(AV)f without the trVptrV^. 

There are only nine terms in Zp that do not have this simplification. They are the 
terms which have at least one power of {AA) on each side, ie. 



0(AA) 



0(AA) 



0(AA) 



0(AA) 



0(AA) 



0(AA) 



0(AA) 



0(AA)' 



+ 



0(AA)' 



0(AA)' 



0(AA) 



0(AA) 



+ 



0(AA) 



0(AA) 0(AA) 
' 0(AAf _j_ 0(AAf 



0(AA) 



D 



D 



0(AA) 



0(AA) 



0(AA)' 



0(AA) 



0(AA) 



0(AA) 



0(AA) 



0(AA) 



0(AA) 



+ 



0(AA) 



0( 



0(AA) 



0(AA)' 



0(AA) 



0(AA)' 



0(AA) 



0(AA) 



0(AA)' 



The evaluation of these terms starts by spliting up the trVptrV^ into Vat V^^ pairs 
around the sides of the plaquette, working out the group integrals on each side (following 
Samuel [^) and then "glueing" the sides together by contracting the spare colour indices 
(ab,cd). 

As an example we give the explicit form for the first of the integrals in ( ]92| ) 



41 
648 
1 

"648 



+ 



324 



^324 

]_ 

162 



_l_ 

'81 



dVidVidVsdVi ymysV^ trVptrV} 
tr{MAi)tr(A2A2)tr(A^A^)tr{MAi) 
[tr(AiAi)tr{A2A2)tr{A^Aj.;A^) + tr(A2A2)tr(A^A^)tr{AiAiAiM) 
+tr(A^A^)tr(AiA^)tr{A^A2A2Ai) + tr(AiAi)tr(AiAi)tr{A2A:iA^A2)] 
[tr{AiAi)tr{A2A^Ai'A^^'A2) + tr{A2A2)tr{A:iAiAi'Ai'A^-i) 
+tr(AzA^)triAiAiA2'A2'Ai'Ai) + tr(AiAi)triAiA2A^'A^'A2'Ai)'^ 
[tr{AiA2'A2'A{)tr{A:iA^iA^) + tr{AiA{A{Ai)tr{A2A^'A^'A2)] 
{tr{AiA2A^Ai'A^^'A2'Ai) + tr{A2A^AiA{A{A^^'A2) 
+tr{A:iAiAiA2'A2'Ai'A^^) + tr{AiAiA2A^'A^'A2'Ai'Ai)] 
tr{AiA2A^Ai)tr(A;A^'A2Ai) . 



(93) 



It it obvious from (|93|) above that the nine integrals are not separable into products of 
link-functions. We re-write Zp = Z° + Z where Z contains the above nine integrals 



minus the corresponding integrals in Z^ (without the trVptrV^). So Jp = 1 + Z'^{Z'p 
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Expanding (Z^) ^ in powers of 0{AA) 



hitriAi^^Ai^^)]^ — ^tr[{Ai^fj_Ai^^Y] + T^[tr(Aj^^ylj_^)] 



1441 

49 
2880 



where we have truncated up to order 0{AA)^ 
Looking at example term in eg. 



1 



0(AA) 



dVidV2dV3dV^—iyif{y2yy3{y4)HrV,trV; or omf 
d4 



0(AA) 



(95) 



0(AA) 



When evaluated this will be an even more complicated function of Ai, Ai, A2, A2, A^, A3, 
Ai, A/^ than (p3D. However although it will have four pairs of '^'^ at sites B and C (so two 
pairs of MM) it will have three pairs of "^"^ at sites A and D which cannot be arranged 
solely into MM pairs. In order to give a non-zero contribution, this diagram will have 
to be multiplied by the factor —^tr^AiA^) coming from (Z^)"^. This will then give a 
diagram with 0{AAY on each side and all \E'\E' fields can be re-arranged into tr{MiMi) 
form. So we must evaluate each of the nine diagrams in (|92D and then multiply each by 
whatever terms in (Z°)~^ (^) which will give a contribution containing pairs of at 
each site. 

There will be only four ways of arranging the 0{AA)'^ terms on each side of the 
plaquette to ensure that on each corner there are either two \E'\E' pairs or four \E'\E' pairs. 
These are shown below. 



0(AA) 




0(AA)' 




0(AA)' 




0(AA)' 


0(AA) 


0(AA) 


or 


o(AA)' or 


0(AA) 1 


0(AA) 


or 


0(AA) 




0(AA)' 




0(AA) 




0(AA)' 



0(AA) 



(96) 



Finally we describe the procedure for summing over all the delta functions which arise 



after substituting < trM^Mf > 



One term which will occur in Zp{Zp 



2Uf5'^^5^^ at each site. 
1 



IS 



0(AA)' 



0(AA)' 



0(AA)" , 



(97) 



0(AA)' 
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where a simple possibility for 0{AA) is tr(Aj^^Aj^^)^. Of course there will be much more 
complicated terms of 0(AA)^ where the A, A are connected by colour indices to other 
A, A on adjacent sides, but the procedure will be the same. 

and this can be rewritten, 

~ (*^^*^"*^^*5^)(7f 7?7f 7?)(*S^*S"*^/*B ), (99) 
or in a more abstract form 

a^.X^^.a^, (100) 

where . represents the contraction over the Dirac indices. So all around the plaquette we 
have [ir(]4iAi)]^[tr(]42A2)]^[*'^(^3^3)]^[i'^(^4^4)]^, written in this abstract form as 

(a^.X^^.a^)(a^.X^^.a^)(a^.X^^.a^)(a^.X^^.a^), (101) 

or 

rr[(a^a^).X^^.(a^a^).X^^.(a^a^).X^^.(a^a^).X^^], (102) 

where Tr is, at the moment, just indicating that all Dirac indices are being summed over, 
(a^a"^) now has 8 Dirac indices, 4 of which it contracts with X^^ and 4 with X^"^. 

The next step, is to rc-arrange the 4 (^'j, ^'j) pairs in (a*a*) into 2 meson pairs. As 
mentioned in Appendix 1, this can be done in three, equally likely ways. The function at 
site B, a^a^ will be of the form 

a^a^ - ^S^^^^^^^^^V^^^J"*^'^^^'", (103) 

where ^5^^^^*^%^^ comes from side AB and from side BC. We 

rewrite this as 

Now we substitute in < Mf'''"^ >= C/jCrf 5°^ to get 

^U^{d^^6^^d''^d^'' + d^''d'^^df^d^^ + d^''6f'^d^^d'^'^] . (105) 

So wc have effectively replaced (a*a*) with ±^UfA'^'^'^^P''^'', with A representing the 
sum of products of delta fns (the bit in {}) , which has 8 Dirac indices. The ± above 
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appears because when we rearrange the we need to be careful about the anti- 

commuting Grassmann numbers. 

Since each Dirac index has 2 possible values (it's a 2x2 representation of the Dirac 
algebra) and half the indices contract to the right and half to the left, we can re-represent 
the function A (a 2^ tensor) with a 16x16 matrix. In (|105|) the four indices I39C,'^ connect 
via four 7 matrices to site A, and pcrxi^ connect via four 7 matrices to site C. We can 
therefore replace each set of four indices (four pairs = sixteen combinations) with a single 
index taking 16 values. So A^^^^t)^'"^^^) = A®^. The same goes for the functions X 



representing the four gamma matrices in ( p.02| ) and we merely need to calculate products 
of 16x16 matrices around the plaquette. The Tr in (|102|) then becomes a matrix trace. 
Actually working out the matrix form for each product of delta functions then becomes 
a matter of book-keeping. 

A helpful result is that, because of the symmetry of the problem and the even number 
of gamma matrices, we can actually replace the Xs with identity matrices. 

We have outlined the procedure for working out one term, where all the corner func- 
tions were of the form (a* a*) but the same goes for each term in Zp{Z'^)~^. In the case 
where there are spare colour indices to contract at the corners we must do this first and 
then re-arrange the ^^!, but each corner will eventually be written like ( tL03| ). 

For the other diagrams in (|96|) there will not be four gamma matrices on the sides 
but rather two (for 0{AA)) or six (for 0{AAY) and these can be rewritten as 4 x 4 or 
64 X 64 matrices respectively. So the corner function will be a 4 x 4 matrix for a corner 
connecting two 0{AA) sides, and for a corner connecting a 0{AA) side and a 0{AAY 
side a 4 X 64 matrix (with the transpose 64 x 4 for an adjacent corner). The general 
procedure however, is identical in each case. Our computation is given in more detail in 
138| and is summarized below. 



1. Work out the group integrals on each side, the result will be a sum of terms like ( |1UU| ) 

2. Write the integral around the whole plaquette as a sum of terms each of the form (|102|) 
contracting all spare colour indices. 

3. Replace each corner functions with the sum of delta functions, being careful to keep 
track of ± due to the anti-commuting nature of (\E'j, ^i) 

4. Rewrite the A functions as 16 x 16, 4 x 4, or 4 x 64 matrices 

5. Trace over the matrices 
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